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C^i . Several approximations are tested by calculating the ground-state energy and the 

' energy of the first excited 0^ state using an exactly solvable model with two symmet- 

, ric levels interacting via a pairing force. They are the BCS approximation (BCS), 

Lipkin - Nogami (LN) method, random-phase approximation (RPA), quasiparticle 
O ■ RPA (QRPA), the renormalized RPA (RRPA), and renormalized QRPA (RQRPA). 

^3 . It is shown that, in the strong-coupling regime, the QRPA which neglects the scatter- 

ing term of the model Hamiltonian offers the best fit to the exact solutions. A recipe 
is proposed using the RRPA and RQRPA in combination with the pairing gap given 
by the LN method. Applying this recipe, it is shown that the normal-superfluid phase 
^ ' transition is avoided, and a reasonably good description for both of the ground-state 

I energy and the energy of the first excited 0^ state is achieved. 

^ ! 

O ■ PACS numbers: 2f.60.Jz, 21.60.-n 

O 
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The random-phase approximation (RPA) [1,2] is one of the most popular method in the theoretical 
microscopic study of nuclear structure. It includes the correlations beyond the mean field models 
^ I such as the Hartree-Fock (HF) approximation with a phenomenological interaction, e.g. the Skyrme 

C3 ' interaction. Being also a computationally simple method, the RPA serves as a powerful tool in 

treating all the excited states in nuclei, which are beyond the reach of the full diagonalization within 
a shell-model basis. The quasiparticle RPA (QRPA), in which the quasiparticle correlations and 
excitations are considered, has shown its prominent role in treating the open-shell nuclei, where the 
superfluid pairing correlations are important. 

Recent developments in nuclear structure including /3-decay and double-beta {(3(3) -decay physics, 
the study of various types of giant resonances, and the prospect of using radioactive beams to study 
nuclei far from the line of /3-stability have sparked off a renewed interest in the efficiency of the RPA 
compared to other microscopic calculations [3]. More refined and exact treatments of the pairing 
problem in nuclei have also been proposed in [4] . 

Recently the accuracy of the RPA in describing the binding energies has been tested using the 
HF-I-RPA calculations within schematic exactly solvable models [5,6] as well as for nuclei throughout 
the sc?-shell and the lower p/-shell [7] . The tests using exactly solvable Lipkin models [5,6] has shown 
that the HF-I-RPA (QRPA) calculations yield a very good ground state energy except for the region of 
the pairing interaction around the point where the pairing gap collapses. Mean while, the tests using 
more realistic shell-model Hamiltonian have pointed out that the binding energy predicted by RPA is 
generally, but not always, satisfactory [7,8]. A number of suggestions have been made to improve the 
reliability of the RPA. They can be classified into two groups. The first group includes the approaches 
to improve the treatment of pairing correlations. Among them are the Hartree-Fock-Bogolyubov 
(HFB) -f QRPA, the number projection [2], the variational method using boson expansion [9], 
and the Lipkin-Nogami (LN) method [10,11]. The LN method is an approximation to remove the 
fluctuations due to the violation of particle number conservation of the BCS approximation, which 
leads to the collapses of the pairing gap at a certain critical value of the pairing interaction. The 
second group includes various treatments of the ground-state correlations beyond RPA, which occur 
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due to the violation of the Pauh principle when the quasiboson approximation is used within the 
RPA or QRPA. In a way similar to the particle number violation within the BCS approximation, 
the violation of the Pauli principle by treating the fcrmion pairs as bosons within the quasiboson 
approximation leads to the collapse of the RPA at a critical value of the interaction parameter. 
Various approaches have been proposed to renormalize the RPA (QRPA) to remove this inconsistency 
[12-17]. 

Since the RPA, first of all, is a theory of excited states, an improvement of the accuracy of the 
RPA should give, as a first step, a better description for both of the energies of the ground state and 
the first excited state simultaneously. In the present paper, we arc going to propose a recipe for such 
improvement using existing approximation schemes, namely the LN method for the pairing gap and 
the renormalized RPA (renormalized QRPA) to treat the ground state correlations beyond RPA. In 
order to compare the result with the exact solution we limit ourselves to the study of the ground 
state and the first excited 0+ state within an well-known exactly solvable two-level model, which 
was introduced for the first time in Rcf. [18], and widely used in the study of pairing correlations. 

The paper is organized as follows. The pairing Hamiltonian applied to the two-level model is 
discussed in Sec. I. The BCS approximation, the LN method, and the results obtained within these 
approximations for the pairing gap in the two-level model under consideration are summarized in 
Sec. II. The RPA, QRPA and their renormalized versions are presented in Sec. Ill, where the results 
for the ground-state energy and energy of the first excited 0+ state are analysed. Also discussed 
in the same section arc results obtained following a recipe which combines the LN-mcthod and the 
renormalized QRPA (renormalized RPA). The last section summarizes the paper, where conclusions 
are drawn. 



II. PAIRING HAMILTONIAN 

We consider the well-known pairing Hamiltonian of nucleons interacting via a pairing force between 
the time-conjugate orbitals with angular- momentum quantum number j 

J jj' 

where ej are the single-particle energies, and G is the pairing constant. The nucleon number operator 
for the j shell is 

= Yl ^]m(^jrn ■ (2) 

m 

The pairing operators Aj and Aj are given by 

4 = y «La^~ , = (4)t , (3) 

V J m>0 

where the tilde denotes the time- reversal operation, e.g. = (— l)-'~'"at__. These operators 



(4) 

-25, y A, . (5) 

The exact solutions of this pairing Hamiltonian have been found in Ref. [19] and known as the 
Richarson's solution. In the present paper, we consider only a simple schematic model. It has N 



satisfy the following commutation relations: 

[A,,A],] = 5,,(l-^^ 
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particles occupying two levels with the same shell degeneracy il = j + 1/2 so that f2 = N/2. The 
upper level simulates the degenerates states with energy e/2 and magnetic quantum numbers m 
(0< m < n), while the lower level is for the states with energy —e/2 and the magnetic quantum 
numbers —m. The distance between two levels is, therefore, e, which will take the value equal to 1 
MeV in all calculations in this paper. The exact solution in this case is easy to be found using the 
SU(2) algebra, which the operators J+ = i/n"^], J_ = (J+)^ and Jo = (Nj - flj)/2 obey. This 
leads to the diagonalization of the Hamiltonian (1), whose matrix elements are 

{J,M[,M^\H\J,M^,M2) = -e{M, - M2)6M„M[SM,,Mi 
-G[2J{J + 1) - Ml (Ml - 1) - M2(M2 - l)]5MuMi5M,,M^ 
-G^J{J+l) - Mi{Mi - l)^/J{J + l) - M2(M2 + 1)5m^ ,M( + l^M2,M^-l (6) 

-G^/J{J + 1) - Ml (Ml + lWJ{J + 1) - M2(M2 - 1)5m, 

with J = fl/2 and — J < Mj < J (i = 1, 2). Among the obtained eigenenergies Si, the lowest 
one, So, is the exact ground-state energy, while the exact energy Uex of the first excited 0+ state is 

Wex = Si — So . 

It is convenient to study a fermion system with superfluid pairing using the Bogolyubov trans- 
formation from particle operators aj^^^ and The quasi- 
particle representation of the Hamiltonian (1) is given in [18,20], which we quote here again for the 
convenience in further discussions: 

H = a + Y, bjJ^j + J2 + A- ) + II dj,, A]Aj. 

3 3 33' 

+ E 33 if) i^l'J^J + ^3-^3' ) + E ^ii' ^-^lA' + -^i'-Ai) + E ^n'^3^3' ■ (7) 
ii' ii' ii' 

Here J\fj is the operator of the quasiparticle number on the j-shell, while and Aj are the creation 
and destruction operators of a pair of time-conjugate quasiparticles: 

^3 = E "U"^- ' A = -7^.Y. '^Im'^]^ > ^3 = iA)^ ■ (8) 
m v 3 rn>0 

Their commutation relations are similar to those for nucleon operators in Eqs. (4) and (5), namely: 

[Aj,A],]=djj,{l-^), (9) 

Wj,A],] = 25jj,Al ' [■^^■"4/] = -2%'^- . (10) 
The coefficients in Eq. (7) are 

a = 2 E ^3^3^j - G(E ^JUjVjf - G E ' (11) 

3 3 3 

i>3 = ^3 ("f - '"'j) + '^■GujVj E ^r'^f'^r + Gv^j ' (12) 

j' 
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Cj = 2^/n'j€jUjVj - G^/TTjiu^j - u|) ^ .J^,UyVy - IG^jUjv] , (13) 

j' 

djj' = -G^JVV (^'4 + «I4 ) = dj'j , (14) 
9j if) = GujVj {u} - 4 ) ' (15) 
G 



hjj, = - vTi/v ("14 + 4^1' ) = ^i'j ' (16) 

fe' = -GujVjUj'Vf = qfj , (17) 

where uj and ?;j are the coefficients of the Bogolyubov transformation. Hereafter we will refer to 
the terms at the right-hand side (RHS) of Eq. (7), which contain the coefficients aj , bj, etc. as the 
a-term, b-term, etc., respectively. 



III. GAP EQUATIONS 
A. BCS approximation 

The BCS equation is usually obtained using the variational procedure to get the minimum of 

the average value of iJ — XN (A is the chemical potential, N = Nj is the particlc-numbcr 
operator) over the BCS ground state taken as the quasiparticle vacuum \0)a, i-c. (Xjm\0)a = 0, 
where |0)q = JJj „>□(% ~^ '^j^lm^''~)\^) with ajm\0) = 0. Within the BCS approximation, only the 

a-term in Eq. (7) contributes, which leads to the well-known BCS equation to determine the gap A 
and chemical potential A: 

A = G ^ n,u,v, , N = Y,^Jl- , (18) 

where the single-particle energy is e'j = ej if the self-energy term —Gvj is neglected, or e'j = ej — Gvj 
if the self-energy term is included. The quasiparticle energy is Ej = ^(e^- — A)^ -|- A^. The Uj and 
Vj coefficients are given as 

The a-term (Eq. (11)) is actually the ground state energy within the BCS approximation, since this 
is the only term that remains in the average over the quasiparticle vacuum \0)a, where the second 
term can be now replaced with —h?/G using the BCS equation (18). 

In the present two-level model, using the BCS equations (18) and the property -\- vf = 1 with 
i = 1 (lower level) or 2 (upper level), it is easy to see that the quasiparticle energy E and the 
chemical potential A are state-independent, namely 

Ei=E2 = E = Gn, ^ = -^- (20) 
The gap A and the u and v coefficients are [6] 
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where 



A = Gf!jl-[^l , (21) 



e neglecting the self-energy term , 

e= { (23) 
2fle/ (2ri — 1) including the self-energy term , 



and 



= ^ . (24) 
The BCS ground state energy in this model becomes 

Ebcs = - 2v'^) - IS?/G - GQ.{u^ + v^) 

-eV2G - (1 - 1 /2Q) /G - GO. neglecting the self-energy term , 

(25) 

-eV(2G(l - 1/20)) - A2(1 - 1/2Q.)/G - GQ self-energy term included . 

This shows that the effect of the self-energy term on the ground-state energy becomes negligible 
at a large particle^ number when 1/2^ <C 1. The major drawback of the BCS method is that the 
BCS wave function is not an eigenstate of the particle number operator N. Therefore the particle- 
number fluctuations AN^ = a{0\N^\0)a - a{\ON\0)l = ^J^j ^ji^jVj)^ = 2A'^/G'^fl [20] cause the 
inaccuracy of this method. Using the exact commutation relation (9), and a(0|A/^ |0)a = 0, we see 
that 

a{0\[Aj,A],]\0)a = , a(0|[^,-,^],]|0)« = , (26) 

which imply that, within the BCS approximation, the quasiparticle-pair operators Aj and Aj behave 
like bosons (the so-called Cooper pairs). Such violation of the Pauli principle between quasiparticles 
has the same origin as the quasiboson approximation used in the RPA, which will be discussed in the 
next section. This leads to the collapse of the pairing gap at the critical value Ga, below which the 
BCS equations (18) yields the imaginary solution. In the present two-level model, Gcr is given by 
Eq. (24). Such kind of critical behavior inspired a speculation of the existence of a phase transition 
from the normal phase, where the gap is absent, to the superfluid phase with a nonzero gap. In 
the absence of the pairing gap, only the sum over the hole (h) state remains in expression for the 
ground-state energy (11), as Vh = 1 and Vp = 0. This gives the HF ground-state energy within the 
present two-level model as 

^HF = (HF|if |HF) = 2 ^ nj.ej, - GJ^ = '^{e + G) , (27) 

jh ih 

where |IIF) = Y\j ^\mW) ^he HF ground state. However, the superfluid-normal phase transition in 
a system with a finite particle number is spurious as it does not exist in the exact calculations [21] as 
well as in the method using particle- number projection [2] , where the gap is finite at all finite values 
of G [2]. Since carrying out the particle- number projection in calculations using realistic spectra is 
numerically complicate, a simple approximate number projection has been proposed, which is known 
as the Lipkin-Nogami (LN) method [10,11] and summarized below for the present two-level model. 
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B. Lipkin-Nogami (LN) method 



The LN method has gained a great popularity as it provides a simple and computationally easy 
way to go beyond the pairing mean-field. Within this method, the particle-number fluctuations are 
removed by adding the term A2-/V^, and carrying out the variational procedure over the average 
value of iJ — XN — A2A^^ in the quasiparticle vacuum \0)a- Details of this method is given in Ref. 
[11]. For the present two- level model, it gives 



(28) 



-' = Ui-t;^), v^-Ul + ^], (29) 



2Gn J 2 V 2Gn 



where 

2Gne 



2Gn + a 



AX2-G . (30) 



Setting the factor a = recovers the BCS equation without the self-energy term, while setting A2 = 
brings back to the BCS equation including this term. The factor a is found substituting (uw)^ in 
the equation for A2: 

4A2 ^ n - 2{uv f 

G 2{2n-i){uvy ■ ^ ' 

This leads to the following cubic equation for a: 

a{2n - l)[{a + 2GVLf - e^] - 2GVLe^ = . (32) 
The ground state energy is given by 

£;ln = ^{u^ - v^) - -pr- - G^[^ - 2{uv f] - 2{a + G)n{uv f . (33) 

Gr 




FIG. 1. Pairing gaps (in units of 2e) as a function of Gfl/{2e) for $1 = 4 (a) and 8 (b). The dotted line 
denotes the BCS gap A obtained including the self-energy term. The dashed line is the BCS result neglecting 
the self-energy term. The solid line stands for the gap Aln given by the LN method. 
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The gaps obtained within the BCS approximation and LN method using Eqs. (21) and (28) for 
= 4 and 8 are plotted against the interaction parameter GQ/(2e) in Fig. 1. With decreasing 
the interaction, the BCS gap A decreases and collapses at Gcr, whose value for the case when the 
self-energy term is neglected is smaller than that obtained including the self-energy term. The LN 
gap Aln, on the contrary, decreases monotonously with decreasing G until G = 0, where Aln 
vanishes, showing no signature of the supcrfluid-normal phase transition. Hence, by removing the 
particle number fluctuations, the LN method erases completely the odd behavior of the BCS gap 
characterized by this phase transition. The system remains in the superfluid phase at all nonzero 
values of the pairing interaction. It is interesting to sec that, due to the suppression of the self-energy 
term —Gv^ by 4A2W^ in the LN method, the BCS gap without the self-energy correction is closer to 
the LN result than the BCS gap including this correction. By comparing Fig. 1 (a) and (b), we also 
see that the difference between two BCS versions (with and without the self-energy term) decreases 
significantly when the particle number increases. 

IV. QRPA (RPA), ITS RENORMALIZATION, AND COMBINATION WITH THE LN 

METHOD 

The RPA ground state includes 2p2h-, ApAh-, 6p6h-, etc. excitations on top of the HF ground 
state. The RPA including the pairing correlations within the quasiparticlc representation is usually 
referred to as the QRPA. The correlations in the QRPA ground state, therefore, are much richer 
than the 2p2h-, ApAh-, etc. correlations in the "diffuse" (quasiparticlc) ground state \0)a created by 
scattering of particle pairs within the BCS approximation. The LN method takes into accoimt some 
2p2/i-corrclations beyond the BCS approximation by using the A2iV^ term. However, this method 
still uses the same BCS ground state |0)q, since it approximately expresses the expectation value of 
the Hamiltonian with respect to the projected state in terms of that with respect to the BCS ground 
state [11]. The discussion below is conducted within the QRPA, where the RPA is obtained as the 
limit when A = 0. Consequently, we summarize also the main features of the renormalized QRPA, 
whose zero-pairing limit is the renormalized RPA. In the present two-level model, the QRPA works 
at G > Gcr, while the RPA is applied at G < Gcr- 

A. QRPA 

The standard QRPA operators, called phonon operators, have the following form in the present 
two-level model 

Qi = Yj"U,) , = (Qt )t , (34) 

j 

The QRPA ground state |0)q is defined as the vacuum for the phonon operator, i.e. Q„\0)q = = 
q{0\QI. The 0+ excited state |0+) is obtained by acting Ql on this vacuum, i.e. |0+) ~ Qt\^)Q- 

The excitation energy u)^ of the state 10+), and the amplitudes X^"^ and vj"^ are found respectively 
as the eigenenergy and the components of the eigenvector of the QRPA equation, which is derived 
from the following equation of motion for the Hamiltonian (7): 

q{0\[SQ,[H,QI]]\0)q=u,q{0\[SQ,QI]\0)q . (35) 

In the standard way of derivation of the QRPA equations, the BCS equation is solved first. Then the 
a- and b-terms in the Hamiltonian (7) is replaced with the BCS result, which is Hbcs = EjMj. 
Using the exact commutation relations (9), we see that, among the remaining terms of (7), which do 
not contribute in the BCS, the d-, h-, and q-terms start to contribute within the QRPA. The c-term 
and g-term do not contribute since, in the commutation with the phonon operators (34), the former 
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gives a number, while the latter leads to the terms ^ AjAj,, ^ AjAj', and A/^ (l— A/'j' /^j'), which are 
left out by linearizing the equation of motion according to (35). Moreover, in order to obtain a set of 
QRPA equations, linear with respect to the Xj'^^ and vj"^ amplitudes, another approximation called 
the quasiboson approximation is made, which implies that the following approximate commutation 
relation holds 

[Aj,A],] = 5jj, , (36) 

instead of Eq. (9). The definition of phonon operators (34) and the quasiboson approximation (36) 
lead to the well-known normalization of the QRPA X^"^ and vj"^ amplitudes 

^[Xj^^xf ) - Yj^^Yj"'^] = 6,,, (37) 

3 

SO that the phonon operators are bosons, i.e.: 

[Q.,Ql,]=5..' (38) 

The quasiboson approximation (36) shows that the quasiparticle-pair operators A^j and Aj behave 
like boson operators when interacting with each other. The effect of Pauli principle represented by 
the last term at the RHS of (9) is just ignored. The set of QRPA equations obtained in this way is 
written in the matrix form as 




The explicit form of the matrices A and B depend on the approximation. Below we compare the 
results obtained within the boson and fermion formalisms. 



1. Boson formalism 
The boson formalism is based on two following assumptions: 

(a) It considers A^j and Aj' as ideal bosons bj and bj, respectively, according to the quasiboson 
approximation (36). The j-shell quasiparticle number operator Mj is then mapped onto a boson 
pair as 

AG=2b]b,. (40) 

This mapping preserves the commutators (10). The Hamiltonian (7) can be then fully expressed in 
terms of the boson operators and hj . 

(b) In deriving the QRPA equations according to (35), the q-term of the Hamiltonian (7) is 
neglected. This term is a special case of the so-called scattering term in the general Hamiltonian 
with a two-body residual interaction. For instance, when the residual interaction is separable, the 
q-term involves a sum of products of two scattering quasiparticle pairs B^.j,{LM)Bj^ji^{LM), 
where Bjj'{LM) = — Y^^^,{j'nni'm'\LM)a^-^a.,~, is the scattering quasiparticle-pair operator. 

The latter is equal to Mj/\/2i + 1 when L = 0. The scattering term is usually omitted in most of 
numerical calculations within QRPA for realistic nuclei in literature. 
The boson mapping of the phonon operator (34) becomes 

Ql E(^i''^bt _ .) . (41) 

3 



8 



The QRPA matrices Ajj' and Bjj' have the simple form: 

^jj' = ^Ej^jj' + djj' , Bjjr = 2hjj' , (42) 
which, in the present two-level model, become 

An=A22 = Gn+^, (43) 



112 = ^21 



2G0 ' 



(44) 



Bi2 = B2i = Gfl-^ . (46) 



2G0 

The QRPA equations give one positive solution equal to 

'^QRPA = 2A , (47) 

while the spurious state associated with the non-conservation of particle number is located exactly 
at zero energy. We see that the energy of the first excited 0+ state above the phonon ground state 
|0)q is just twice the pairing gap, the same for the lowest two-quasiparticle excitation above the 
quasiparticlc ground state \0)a of a system with an even particle number. 

In the normal- fluid phase (A=0), the quasiparticlc operator ajm becomes the particle- (p) creation 
or ho\e-{h) destruction operator depending on whether the level is located above or below the Fermi 
level, namely 

4m - ( !! ^ ^ (48) 

Therefore, the following boson mapping for operators and Aj holds 
-t ^ / b] if > A r b, if ej > X 

^^^{hj ife-<A ' lb] ife,<A ' ^^^^ 

The boson mapping for the number operator Nj is given as 

\2(n,-btb,) ife,<A ' 

which preserves the commutation relations (5) and the particle number on the j-shell equal to 20, j. 
The RPA matrices for the two-level model are 

An = A22 = 2(e - GO) , A12 = A21 = . (51) 

Bn = B22 = , B12 = B21 = -2Gn . (52) 
They lead to the RPA phonon energy 
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(53) 



The critical value of Ga where the RPA collapses is the same collapsing point of the BCS obtained 
without the self-energy term (Eqs. (23) and (24)). Therefore, in order to have the superfluid regime 
start at the same critical point, the gap A in Eqs. (43) - (47) should be calculated neglecting the 
self-energy term in the BCS equation (21). This gives 

The phonon energies (47) and (53) are exactly those obtained for the first time in Ref. [18] using 
the space- variable technique and the gap A^^^ (54). 



2. Fermion formalism 

The fermion formalism does not use directly the quasiboson approximation (36). Instead, it 
employs the exact commutation relations (9) and (10) to rearrange the results of calculating the 
commutators [-ff, ^j] and [if, into the normal order. Then, in the process of linearizing the 
equation of motion (35) the following "average" quasiboson approximation is used 

q(0|[A-,4]|0)q = '^^-^' ' (55) 
assuming that the quasiparticle occupation number tij in the RPA ground state is zero, i.e. 

^ _q(0|AG-|0)q 



2% 



= . (56) 



The QRPA matrices are obtained in this way as 

n 



Ajf = 2{Ej + 2qjj)6jj, + djj, , Bjj, = 2(1 - ^5jj. )hjj, . (57) 

3 



Their explicit form in the two-level model is given as 



112 = ^21 



2Gn ' 



(59) 



A2 A2 

^"=^^^ = "2Gn^ + 2GO ' (^°) 

Bi2 = B.=Gn-^. (61) 

The BCS equation (21) for the gap A includes the self-energy term. The 4gjj(5jj/-term in the 
expression for Ajj (Eqs. (57) and (58)) appears due to the use of the exact commutation relation 

(10) to calculate the commutator between the q-term of the Hamiltonian (7) and ^] (or Aj) as 
follows: 
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E Wi.^ii > -4]] = 4fe-^] + 4 E dn'A^i' ■ (62) 

jij'l i' 

The last term at the RHS of (62) does not contribute in the hnearization within the QRPA. The first 
term leads to the above-mentioned Aqjjdjji-term. It is worthwhile to notice that, if one rearranges 
the quasiparticle operators in the q-term of (7) to the normal order as 

jj' j jmj'm' 

and then drops the last term at the RHS (although there is no rigorous justification for doing so), 
the remaining term in the commutation with gives 

E i^^^^r. ' A'] ^ E i-^o' , ^J] = ^UJ^] , (64) 
jlj'l f 

i.e. just the half of the first term at the RHS of Eq. (62). This explains the difference between Eq. 
(58) above and Eq. (38) of Ref. [6] by a factor of 2 in the denominator of the second term at their 
RHS. 

The energy of the first 0+ excited state w obtained solving the QRPA equation with the matrices 
(58) - (61) has the form 



The energy of the spurious state in this case is imaginary. If one uses the approximation (64) instead 
of (62), the positive solution is given by Ref. [6] as: 



while the energy of the spurious state is exactly zero. Neglecting the q-term in (7) yields the energy 
of the first 0+ state as 



^o = 2AJ(l-l-](l + -^] . (67) 



4oy V 4G2r23 

There are two features of the fermion formalism, which are worth noticing here. The fact that in 
the process of linearizing the equation of motion, the exact commutation relations (9) and (10) have 
been used to calculate the commutators [H, A^j] and [H, Aj] does not reduce matrices (57) by setting 
Qjj' = to those given within the boson formalism (42). The difference in Bjj still remains. Another 
feature is that, when the q-term in (7) is omitted, the spurious mode in the fermion formalism is 
shifted to a positive value equal to 



However it is much smaller than the pairing gap A especially at large il, therefore, remains well- 
isolated. Comparing Eqs. (65) - (67) with the boson energy Wqr^pa (Eq. (47)), it is easy to see 
that 

(Ji < woi < Wo < i^QRPA — 2^ • (69) 
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The ground-state energy is now calculated following Refs. [1,6] as 



Sqrpa = Sbcs + (^11 + ^22)] . (70) 

In the limit of zero gap, the pp-KPA equation is obtained from the QRPA equation discussed 
above putting A = 0, v = vi = U2 =1, u = ui = V2 = 0- The solution of this RPA equation is 
decoupled into the addition and removal modes, which have been discussed in detail in Ref. [2,6]. 
In the present two-level model, these two sets of equations can be written in one matrix equation as 



All Bi2\fRp'\ fl 0\ fR]}^^ 



LOr, (r = 1, 2 ) , (71) 



where 



An=€ + 2\ + 2G-Gn , ^22 = e - 2A - Gf2 , B12 = Gfl , (72) 
Wr are found as 



and = Xa , R^f^^ = , R^^ = Y,- , Rf^^ = , wi = Wa , and (J2 = — Wr • The energies Wa and 



Wa = -G - 2A ± ^/{e + G){e + G-2GQ) . (73) 



uJr = G + 2X± V(e + G)(e + G-2G0) , (74) 

where only the sign + in front of the square roots corresponds to the positive values for these 
energies. Here we still keep the factor ±2A, which is useful to connect the RPA solutions with the 
QRPA one at the critical point G = Gcr, where 2 A = — Gcr according to Eq. (20). 

In order make a comparison with the boson formalism, where there is only one boson state, we 
apply the sum-rule method, representing the RPA phonon operator as 

Q^ = Qi + Ql. (75) 

The RPA ground-state wave function |RPA) can be written as a direct product of the ground-state 
wave functions of the (orthogonal) additional and removal modes 

|RPA) = |0)a ^ \0)r ■ (76) 

Using the Thouless theorem for the energy-weighted sum rule <Si [2] with respect to the phonon 
operator (75) 

<5i = ^4'Ll(^il^3^|RPA)|' = ^(HF|[gt, [if,Qt]]|HF) , (77) 

i 

it is easy to sec that the LHS of Eq. (77) is equal wrpa since = Q^'jRPA) {i = 1), while the 
RHS is equal to 



(HF|[Qt, [if,Qt]]|HF) + (HF|[Qt, [H,Ql]]\RF) 



Wa|(!^a|Q];|RPA)|^+a;r|(!^r|Qj|RPA)|' =a;a|(!^aka)|'+Wr|(!/r|l^r)|' =a;a + Wr . (78) 

(The crossing terms (HF| [Q} , [H, Q,l]] |HF) {i i') vanish as can be easily checked using = Sw). 
Therefore 
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WrpA = + t^r • (79) 



The ground-state energy is then given by 



Erpa = ^^HF + 2 (WRPA - ^11 - A22) = V(e + G)(e + G-2Gr!) - (e - + G) , (80) 

which is exactly the same expression obtained previously in Eq. (45) of Ref. [6]. 

The energies of the ground-state and the first 0+ state obtained within the fermion formalism 
of RPA and QRPA for fl — A are compared with the exact energies in Fig. 2. The figure shows 
that, in the superfluid regime, except for the region close to the critical point G — Get where the 
QRPA collapses, the approximation (64) fits the exact result for the ground-state energy better 
than (62). For the energy of the first 0+ state, both of the QRPA versions, which include the q- 
term of the Hamiltonian (7) give practically the same values. However, they are obviously smaller 
compared to the exact solution. This discrepancy increases with increasing the interaction. From 
Eqs. (65) - (67) we found that, in the limit of infinite G the solution cjq (q-term neglected) becomes 
2Ay'l — l/4r2 ~ 3.88A, while the solutions uji based on the approximation (62), and wio based on 
(64) [6] are equal to 0.96wo and O.QSwqj respectively. The only QRPA approximation that fits well 
the exact results for both of the ground-state and excited-state energies the one, which neglects the 
q-term of the Hamiltonian (7) . The results within this approximation practically coincide with the 
exact ones at large G. The RPA energy wrpa — 00^ + of the first 0+ state exhibits the well-known 
behavior. It decreases with increasing G and collapses at the same critical point G = Gcr, from 
which the the normal-fiuid phase ceases to exist, and the superfiuid phase begins. Meanwhile, the 
exact solution for the first excited state has only a bending in this region, showing no signature 
of such phase transition. For the ground-state energy the exact result shows a completely smooth 
curve, while the critical point is clearly seen in the approximations. 




FIG. 2. The ground-state energy (a) and energy of the first excited state (b) (in units of 2e) as a function 
of G'f2/(2e). The thick solid line is the exact result. The thin solid line shows the QRPA result neglecting the 
q-term in the Hamiltonian (7). The dashed line denotes the QRPA result of Ref. [6]. The dotted line stands 
for the result obtained using the matrices (58) - (61). The dash-dotted line represents the RPA result. 

As shown in Fig. 3 for = 4, the QRPA energy obtained without the q-term of (7), and 
the RPA energy wrpa also match well the solutions of the boson formalism, especially after shifting 
Gcr = Gct"* in the latter to the value of e/ (257 — 1) used in the fermion formalism due to the inclusion 
of the self-energy term (Eq. (24)). 
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FIG. 3. The energy of the first 0"*" state (in units of 2e) given by the boson formalism (dashed hne) in 
comparison with those obtained within the QRPA (thin solid line) neglecting the q-term of the Hamiltonian 
(7), RPA as the sum + (dash-dotted line), and the exact solution (thick solid line) for = 4. In (a) 
the boson energy [Eqs. (47) and (53)] has been obtained using Gcr = = e/2Q (neglecting the self-energy 
term). In (b) the boson energy has been calculated using the value Gcr = e/(2f2 — 1) obtained including the 
self-energy term. 



B. Renormalized QRPA 



The collapse of the BCS approximation and QRPA (RPA) has the same origin of neglecting 
the Pauh principle between quasiparticle pairs operators in the BCS approximation (26) and the 
quasiboson approximation strictly (36) or in average (55). The LN method approximately corrects 
this inconsistency within the BCS approximation. For the QRPA this is done by the renormalized 
QRPA. 

The essence of the renormalized QRPA is to replace the quasiboson approximation in the form of 
Eqs. (36) or (55) with the average value of the commutator 

(0| [A, , ^],] |0) = D,S,,, , D,^l^ 2n, , (81) 



in a new ground state |0), where the correlations beyond the QRPA due to the fermion structure of 
\] and Aj 



the quasiparticle pairs J^: and Aj are taken into account, namely 



:^miO)^0, (82) 



2f2j 



instead of the assumption (56). The phonon operators are renormalized as 

S. = E -^i^^^^A - ^A) . 2. = iQl)^ . (83) 



so that the condition for phonons to be bosons within the correlated ground state |0) 

(0|[Q.,Q^,]|0)=<5.., (84) 

leads to the same normalization condition for the amplitudes x'^-'^^ and y'"^^ as that of the QRPA, 

i-^- '^ji'^j'^^'^j" ^ ~ •^j'^^-^j'^ ^) ^ ■ '^^'^ factor Dj is calculated according to the approximation 
in Ref. [16] as 
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The renormalized-QRPA matrices Ajji and Bjj' arc given as 

Ajf = ^Ej + 2qjj)5jj, ^f'Qj'j" (1 - Dj„) + Djdjj, , (86) 



Bjj>=2{D,--^^djf)hjj, . (87) 

Since only the omission of the q-term in (7) within the QRPA reproduces well the exact solution in 
the present two-level model at large values of G, we discuss below only the case when qjj' =0. The 
renormalized QRPA matrices in this case read 

An = A22 = G17(2-D) + ^ , = ^21 = -^(4 - 32?) , (88) 

The energy of the first excited 0+ state is found as 
'^RQRPA = 2A4 



40 y A2 



By setting Dj = 1 {nj = 0) in Eqs. (88) - (90), the QRPA limit in Eqs. (58) ~ (61) (without 
the second term at the RHS of (58) since the q-tcrm is neglected), and ujq given by Eq. (67) are 
recovered . The ground-state energy is calculated using Eq. (70), using the energy wrqrpa given 
by Eq. (90) and the matrices An = A22 given by Eq. (88). 
The renormalized RPA matrices {G < Gct) are given as 

Au=e + 2G- DGQ. , A22 = e - DGQ. , A12 = A21 = , (91) 

Bn = B22 = , B12 = B21 = DGQ. . (92) 
The positive renormalized RPA phonon energies are found as 



w^ = -G-2\+ ^{e + G){e + G-2DGQ) . (93) 



u, = G + 2\+ ^/{e + G){e + G-2DGn) , (94) 

WRRPA = + = 2V(e + C?)(e + G-2£)GO) , (95) 
and the ground-state energy is 

£;rr,pa = y/{e + G){e + G~2DGQ) - (e + G - DGQ) . (96) 

As has been shown in Ref. [16], the presence of the factor D in the renormalized RPA matrices 
reduces the interaction in such a way that the critical point where the RPA collapses is completely 
washed out. The energy 2a, uJr, and consequently, tiiRRPA are always real. As for the renormalized 
QRPA, the presence of the pairing gap A makes it still collapse if A is calculated within the BCS 
(Eq. (21), but it is no longer the case if the LN pairing gap Aln (28) is used. 



15 



C. LN- method + renormalized RPA (renormalized QRPA) 



We have seen in the preceding sections that the LN-method aUows us to avoid the phase transition 
of the BCS approximation, the QRPA without the q-term in (7) provides us with the best fit of the 
exact resuhs for the energies of the ground state and the first excited 0+ state, while the renormahzed 
RPA is known to smear out the coUapse of the RPA due to the violation of the Pauli principle within 
the quasiboson approximation [16]. Therefore, in order to avoid the phase transition and to give 
at the same time a good description for both of the energies of the ground state as well as the fist 
excited 0^ state, we propose here the following recipe. 




FIG. 4. The energies of the ground state (left panels) and of the first 0"*" excited state (right panels) (in 
units of 2e) at several values of Q. The thick solid line is the exact result. The dotted line is the QRPA 
result using the matrices (58) - (61). The dash-dotted line denotes the RPA result. In (a) - (c): the dashed 
line shows the LN result, while the thin solid line represents the ground-state energy calculated according to 
(ii) of Section IV C. In (d) - (f): the dashed line stands for renormalized QRPA result using the LN pairing 
gap and neglecting the q-term in (7); the double-dash-dotted line represents the renormalized RPA results; 
the thin solid line is the energy ljq+ according to Eq. (97). 

i) The QRPA equations (without the q-term in (7)) are solved using the pairing gap Aln (28) 
found in the LN method. 
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ii) The ground state energy is calculated using Eq. (70), in which the LN pairing gap Aln is used 
instead of the BCS gap A (21) and the renormalized QRPA matrices given by Eq. (88) are used 
instead of Eq. (58). 

iii) The first excited 0+ state is presented as a mixed state of those obtained within the renor- 
malized RPA and renormalized QRPA (in (i)). Applying the sum-rule method and following a 
derivation similar to Eqs. (76) (78), wc find that the energy ujq+ of this mixed state is the sum of 
the renormalized RPA and renormalized QRPA energies, namely 

UJo+ = Wrrpa + Wrqrpa , (97) 

where cjRRPA and ojqrpa arc given by Eqs. (95) and (90), respectively, with Aln used in the latter. 
The results obtained using this recipe will be referred to as the combined results. 

The combined results for ground state energy and the energy of the first excited 0+ state are 
presented as thin solid lines in Fig. 4 at several values of in comparison with the exact results and 
those of RPA, QRPA and their renormalized versions. The result obtained within the renormalized 
RPA and the combined result (i) for the energy of the first excited O"*" state show that the superfluid- 
normal phase transition is completely washed out. For the ground-state energy, the combined result 
(ii) is closer to the exact one compared to the LN result in the weak-coupling region (G <C Gcr) 
and in the region around Gcr- It coincides with the LN and exact results in the strong-coupling 
region (G :» Gcr) (See the thin lines in Fig. 4 (a) - (c)). The combined result (iii) agrees reasonable 
well with the exact results (See the thin lines in Fig. 4 (d) - (f)), especially for larger O. Hence, 
this numerical test confirms the validity of the assumption (iii) that the exact excited 0^ state 
may be considered as a mixed state where both of the normal and superfluid phases coexist at 
G ^ 0. The weak-coupling region is dominated by the normal-fluid phase, where the energy u)o+ 
of the first excited 0+ state decreases monotonously with increasing the pairing interaction G. The 
strong-coupling region is dominated by the superfluid phase, in which ujq+ increases with increasing 
G. 



V. CONCLUSIONS 

In this work, a schematic model of two symmetric levels interacting via a pairing force has been 
used to test several well-known approximations by calculating the ground-state energy and the 
energy wo+ of the flrst excited 0+ state. Results obtained within the BCS approximation, LN 
method, boson and fermion formalisms for RPA and QRPA, renormalized RPA, and renormalized 
QRPA have been analyzed in details. It is shown that the common version of the QRPA, which 
neglects the scattering term (the q-term of the Hamiltonian (7)) and the boson formalism give the 
closest results to the exact ones for both of the ground-state energy and u)(,+ at the values of the 
pairing interaction parameter G ^ Gcr- Meanwhile the QRPA version used in Ref. [6] as well as 
the QRPA using the exact commutation relations between quasiparticle-pair operators to treat the 
scattering term do not describe well the energy of the flrst excited 0+ state at large G. 

A recipe has been proposed which combines the results obtained within the LN method for the 
pairing gap Aln, the renormalized RPA, and the renormalized QRPA neglecting the last term (q- 
term) of the Hamiltonian (7). The combined results agree reasonably well with the exact ones for 
both of the ground-state energy and Wo+, showing no signature of a sharp superfluid- normal phase 
transition. The agreement is better at a larger particle number. The results suggest that the exact 
excited 0^ state can be decomposed approximately into two components, which correspond to the 
normal and superfluid regimes, respectively. The weak-coupling region is dominated by the normal- 
fluid phase, while the strong-coupling region - by the superfluid phase. Since the proposed scheme is 
based on rather simple, well-known, and numerically accessible approximations, its future extension 
toward an application in realistic nuclei may be useful. 

For a more self-consistent approach, one can derive the set of RPA equations using the Hamiltonian 
(7) with the coefiicients Uj, Vj, X^"^ and vj'^^ left as variational parameters when minimizing the 
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average energies over the ground state and the first excited 0+ state. This may serve as a goal of 
the future study. 
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